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Abstract 

Using a parameterisation of general self-adjoint boundary condi- 
tions in terms of Lagrange planes we propose a scheme for factorising 
the matrix Schrodinger operator and hence construct a Darboux trans- 
formation an interesting feature of which is that the matrix potential 
and boundary conditions are altered under the transformation. 
We present a solution of the inverse problem in the case of general 
boundary conditions using a Marchenko equation and discusss the spe- 
cialisation to the case of graph with trivial compact part, ie. diagonal 
matrix potential. 



1 The matrix Schrodinger operator on the semi- 
axis 

We consider here the matrix Schrodinger operator on the semi-axis, ie. 

£^-^ + Q(x) 

on L^(M+, C") where = [0, oo). The potential matrix Q{x) is assumed to 
be hermitian, absolutely integrable with absolutely integrable first moment 
and continuous on the open semi-axis [IJ. The matrix Schrodinger operator, 
defined on smooth functions vanishing at the origin and with compact 
support, is a symmetric operator with deficiency indices (n, n). Using von 
Neumann extension theory [2] we may parameterise all self-adjoint exten- 
sions of this operator by unitary mappings between the deficiency subspaces, 
ie. U G U(n). For practical purposes, however, it is more convenient to de- 
scribe the self-adjoint extensions in terms of self-adjoint boundary conditions 
at the origin. 

It is well known that the construction of self-adjoint extensions is analogous 
to the description of Lagrange planes in a hermitian symplectic space [3] . In 
[3] these Lagrange planes are parameterised in terms of two n x n matrices 
such that their product is hermitian. Really, in the case of a hermitian sym- 
plectic space which admits a canonical basi£], the Lagrange Grassmannian 



^By our definition, a hermitian symplectic space for a symmetric operator with non- 
equal deficiency indices does not admit a canonical basis or Lagrange planes [3]. 
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is isomorphic to the unitary group U(n) and we are able to expHcitly pa- 
rameterise the Lagrange planes, and hence self-adjoint boundary conditions, 
in terms of a unitary matrix — for details see jSj. In the case of the matrix 
Schrodinger operator on the semi-axis the self-adjoint boundary conditions 
at the origin are given by 



(1) 



Then the solution of the matrix Schrodinger equation £H = AS with bound- 
ary values 

1 1 

(2) 



satisfies these boundary conditions. Using the Jost solutions, F±, solutions 
of the homogeneous equation CF± = XF±, with asymptotic behaviour 

lim F±{x,k) ^ k = VX 

x— >oo 

and no prescribed behaviour at the origin, we can write 

E{x,k) = F_{x,k)M_{k) + F+{x,k)M^{k). (3) 
We define the scattering wave solution 

^'(x. A;) = E{x, k)MZ^ = F_ + F+S{k) 

where S{h) is known as the scattering matrix. The coefficients M_|_ can be 
evaluated by taking the Wronskian of H and or F_ [3] 



M. 



± 



'2ik 



(4) 



where F-|-(A;) = F±{0, k) are known as the Jost functions and ^ is the invo- 
lution Y\x, k) = Y*{x, k). The Wronskian of and H 



Wi-ElE] 



A*B - B*A = 0, 



is always zero. Moreover, if we write H in terms of the scattering wave 
solutions 

W{r), H} = m1t^{f1 + S^fI,F. + F+S}M_ 



2ikM^ 



M 



0. 



we see, since = S* for A; € M, that the scattering matrix is unitary for 
real k. 

If we diagonalise U, and use the well known asympototics of the Jost func- 
tions [H [3] in the above expression for , we see that the scattering matrix 
has the following asymptotic behaviour: 
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Lemma 1.1 Given the self-adjoint operator L, with associated unitary ma- 
trix U defining the boundary conditions of C, the scattering matrix of C has 
the asymptotics 

lim S{k) ~ U 

fc— >oo 

where U is a unitary hermitian matrix Lf = U* derived from U by applying 
the map 

1 :zGT\{-1} 
^ -1 :z = -l 

to the spectrum of U . 

Here T is the unit circle in C. The matrix since it is hermitian, defines 
projections 

which may be used to define a factorisation of the original operator C 

Theorem 1.2 Given the self-adjoint operator C we can formally factorise 
it as 

jC = D*D 



where 



D 



dx 



and the functions in the domain of D satisfy the following boundary condi- 
tions at the origin, 

Vio = 0. 

Furthermore, V is a hermitian matrix which satisfies the Riccati equation 
with the potential Q{x) on the right hand side and has initial value satisfying 



PV\o 



-PH 



(5) 



where H is a bounded hermitian matrix specified by the boundary conditions 
ofC. 

Proof 1 The fact that V{x) satisfies the Riccati equation is well known. We 
diagonalise the matrix U. Suppose that I of the eigenvalues of U are —1, 
then we choose the first I elements of the basis used in the diagonalisation to 
be the eigenvectors with eigenvalue —1. In this basis the boundary conditions 
of C can be written 



-il 







-tan((^/2) 



V'lo + 





I(n-0 



x\0 







where ip is the boundary value (under the change of basis) and (p is the 
(n — I) X (n — I) diagonal matrix with entries which are the eigenvalues ofU 
excluding the eigenvalues —1. This means that we can write the boundary 
conditions for jC as 



P^P^ + PHi^l^ = 



(6) 
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where H is the hermitian matrix depending on the original boundary condi- 
tions. Consider the operator D*D. From D we get the boundary condition 



and from D* we get the boundary condition 

P[Di;]\, = 

which is 

PVxlo-^^loV'lo = 0. 

So we see that as long as the initial value of V satisfies equation D*D 
has the required boundary conditions. 

It is well known that the Riccati equation can be linearised to the Schrodinger 
equation at zero energy, let us denote by Hq these zero-energy solutions. 
Then it is natural to ask how we can express the coefficient V{x) in terms 
of Hq. The following theorem is proved in [3], we merely quote it here. 

Theorem 1.3 The hermitian matrix V can be written 

V{x) = Eo^x{x)Eq^{x) 

where Ho(x) is the matrix of solutions of the Schrodinger equation at zero 
energy satisfying the boundary values of equation specified by the unitary 
matrix Uq and subject to: 

1. The matrix Uq satisfies 

PUo = PU. (7) 

2. The potential Q{x) is continuous in some neighbourhood containing 
the origin. 

The inverse problem for the matrix Schrodinger operator on the semi-axis 
with Dirichlet boundary conditions (ie. vanishing of the functions at the 
origin) is described in Agranovich and Marchenko [ij. In [3j the author ex- 
tends Agranovich and Marchenko's result to arbitrary self-adjoint boundary 
conditions at the origin; the inverse problem may be reduced to a Marchenko 
equatiorll 

POO 

G{x + y) + K{x,y)+ / K{x,t)G{t + y)dt = Q x<y, (8) 

Jx 

where 

G{t) = y Cfe-^'' + — / {S{k) - U)e'^'dk. 

27r /„no 

Here \i = —k[ are the discrete eigenvalues and the Ci are non-negative her- 
mitian matrices known as the normalisation matrices. The set {^(A;); Ki,Ci,l = 
1, . . . , A^} is known as the scattering data. As this is a lengthy derivation 



^It is also possible to reformulate it as a Riemann-Hilbert problem [3]. 
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we will not repeat it here. If we are able to solve equation ([8]) for the kernel 
of the transformation operator K{x, y) we can recover the potential matrix 
from the well known identity QJ 

We can also recover the self-adjoint boundary conditions at the origin from 
the inverse problem via 

which follows from the definition of the scattering wave solution plus equa- 
tion ([2]). Consequently, the solution of the inverse problem allows us to 
recover not only the potential but also the self-adjoint boundary conditions 
at the origin. 



2 The Schrodinger operator on the graph with 
trivial compact part 

The motivation for studying the matrix Schrodinger operator is that, in the 
case of diagonal potential, it may be identified with the Schrodinger opera- 
tor on the non-compact graph with trivial compact part — here we mean the 
graph consisting of n semi-infinite rays with the origin of each ray identified 
with the single vertex of the graph. Although these are really two different 
operators, for the purposes of the inverse problem they may be identified: 
each component of the vector fuction on which the matrix Schrodinger op- 
erator acts is identified with the value of the function on one of the rays of 
the graph. 

In this case the self-adjoint boundary conditions at the origin play a crucial 
role; for instance if we have Neumann or DirichlelH boundary conditions 
there is no interaction between the rays and the graph decomposes into n 
semi-axes for which the solution of the inverse problem is well known. The 
self-adjoint boundary conditions at the origin describe the interaction be- 
tween the rays and they may also be thought of as inducing a 'zero-range' 
potential at the origin [3]. It is for these reasons that we consider the matrix 
problem with general boundary conditions above. 

For the problem on the graph — the diagonal matrix potential — the matrix of 
Jost solutions is clearly a diagonal matrix and so too is the kernel K{x,y). 
As a result the Marchenko equation ([8]) has along the diagonal n scalar, 
independent (in the sense that each diagonal entry of K{x, y) appears only 
once) Marchenko equations. Consequently, the inverse problem on the graph 
decouples to n scalar inverse problems which can always be solved using 
only the n diagonal entries of the scattering dat^. Using this scattering 
data we can recover the potential on the rays plus the self-adjoint boundary 

^Agranovich and Marchenko fT| consider the inverse problem for the matrix operator 
on the semi-axis, however, they only consider Dirichlet boundary conditions which is not 
interesting in the case of diagonal potential. 

* Actually, from equation ^ we need the diagonal entries of S{k) and Cf . 



5 



conditions/zero-range potential at the origin. 

If the self-adoint boundary conditions at the origin are specified it may be 
possible to recover the potential on the rays using a smaller set of scattering 
data. For Dirichlet or Neumann boundary conditions there is no interaction 
between the rays and we need all n elements of the scattering data in order 
to recover the potential. On the other hand, if there is an even number of 
rays n = 2m and we have boundary conditions at the origin so that the 
graph decomposes into m copies of the whole real axis it is well known that 
we can recover the potential using only m reflection coefficients and normal- 
isation constants [5J. 

Let us consider flux-conserved boundary conditions which are defined by 
continuity at the origin 



>1|0 — <P2|o — • • • — (Pn\o , 



E'^1 = (9) 



i=l 



plus conservation of flux. These are self-adjoint and we provide a brief proof 
that in this case we need only n — 1 of the diagonal elements of the scattering 
data in order to recover the potential (for details see [3]). We are able to 
show that the dispersion function is equal to [3] 

jn-l " f . 

det(2iA;M_ (k)) = M{k) = f+,if+,2 • • • f+,n V 7^ 

where f±^j denotes the Jost function for the z-th ray of the graph and 
f+ = liiUx^o df+^i{x, k)/dx. The scattering matrix has entries [3j, 

Suppose we are given the first n — 1 diagonal entries of the scattering ma- 
trix for real k, the reflection coefficients, Rj{k) = Sjj{k), plus the discrete 
eigenvalues ki = ini and the first n — 1 diagonal entries of the squares 
of the normalisation matrices, denoted hij, for each eigenvalue. Since the 
Marchenko equation degenerates into n independent scalar equations on the 
rays we can use this data to recover the potential and Jost solutions on the 
first n — 1 rays. Equation (jlOj) implies 



Sij 



for i 7^ i. 



Consequently, we have enough information to recover the (n — 1) x (n — 
1) minor of the scattering matrix formed by deleting the last column and 
row. Since S{k) is unitary for /c G M we can recover the magnitudes of the 
remaining entries of the scattering matrix. Let us consider any entry in the 
last column not on the diagonal 



Si', 



2i«A;f+,if+,2---f+,„-i 



nf+,jM(fc) 



i ^ n. 
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As we have the Jost functions for i = l,...,n — Iwe can solve this for the 
magnitude of the dispersion function on the real axis. Now M{k) is analytic 
in the upper half-plane [3j, has known magnitude on the real axis and known 
zeroes — the discrete eigenvalues — and so we can recover M{k) in the upper 
half-plane. To do this we consider the 'normalised' dispersion function 

1 ^ 

Here iki are the zeroes and mi the orders of the zeroes of M{k). In [3] we 
describe in detail how, using the bi^i, to find the orders mi and also show 
there that, subject to the absence of vitual levels, A4{k) is bounded and 
non-zero on the closed real axis with asymptotic 

lim M{k) = 1. 

|fc|-^oo 

Furthermore, by definition M{k) is analytic in the upper half-plane with no 
zeroes there. By a simple application of the Cauchy integral formula and 
the Plemelje formula (see [3] for a similar calculation in the scalar case) 

argM(A;) = 



where we take the principle value of the integral. Consequently we recover 
M{k). From M{k) we can use equation (jlOp to recover Sm for i ^ n. We 
now have enough information to recover the potential on the last ray: from 
the scattering matrix we can find f+^n and from M{k) we can find f+_n- 
These two functions are enough to recover the potential on the ray as they 
provide the scattering data for the Schrodinger operator on the semi-axis [T] 
(with Dirichlet boundary conditions in the cited text). 
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k + iKi 
k — iKi 



N 



lni"(A; + i) + mi In 



1=1 

ln[\M{k')\/\k' + i\] 
k' -k 



k — iKi 
k + iKi 

dk' 
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